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—j , Possibilities for detecting heavy Majorana neutrinos (iV's) at future e/j, collid- 

^sO ', ers are investigated. In contrast to the e~e + colliders (LEP200 and NLC), the 

O ' 



center-of-mass (CMS) energies achieved at e/x colliders can be much higher 
and the Z-mediated s channel is excluded automatically. This opens the at- 
c| ■ tractive possibility of having high production cross sections for iV's and at 

Mh, the same time probing only the strength of charged current couplings of iV's 

Qh! (NWe and NWpi). The production cross sections and the expected num- 

bers of events for the reaction e^/x — > NN — ► W i^W f^ are calculated 
for various masses M of the Majorana neutrinos and for the CMS energies 
y/s = 0.5-6.0 TeV. The values of the charged current coupling parameters are 
set equal to their present upper bounds. We obtain reasonably high produc- 



tion cross sections. Further, the effects of the off-shell intermediate iV's turn 

out to be significant only at y/s ~ 2 TeV. 
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I. INTRODUCTION 

One of the basic questions in high energy physics is: Are neutrinos Dirac or Majorana 
particles? In the absence of right-handed currents, it is virtually impossible to discern 
the nature of the light neutrinos [p]]. However, if heavy neutrinos (M ~ 10 2 GeV) exist, 
then present and, even more so, future experiments could establish whether such neutrinos 
are Majorana or Dirac. Many theoretical investigations have been carried out |2| into the 
production cross sections of heavy Majorana neutrinos (iV's) as predicted by various models, 
via e~e + , ep and pp collisions, and mostly within specific models or classes of models. As a 
rule, it has been assumed that the center-of-mass (CMS) energy y/s in the process is high 
enough for the production of on-shell (OS) heavy iV's (e + e~ ,ep,pp — > NN). Further, the 
effects of the off-shell (nOS) iV's have been ignored. We note that the detection of NN 
pair production events can be realized by identifying the final decay products (on-shell) 
of the Majorana neutrinos - primarily the leptons £/s and the W's. Thus, the reactions 
£-£'+ -> NN -> W ± £fW ± £J -> jets + £f£f, which violate the total lepton number, would 
be a clear signal of the Majorana character of the intermediate neutrinos (the latter need 
not be on-shell). Here, £ and £' are the two initial leptons, either e~e + or e/i. 

In a previous paper || , we demonstrated the importance of taking into account the con- 
tributions of off-shell intermediate iV's in such processes, in e~e + colliders such as LEP200 
and future Linear Colliders. Several advantages of efi colliders were pointed out by the 
authors of ||. In the present paper, we present some results of calculations for the afore- 
mentioned reactions for efi colliders. The main advantages of the efi colliders, in contrast 
to the e~e + colliders, from our perspective, are the following: 

• /i's can easily be accelerated to very high energies due to reduced synchrotron radiation 
loss (since m^^>m e ); 

• the s channel (Z-mediated) is not present now. 

Due to the first point, the event rates can be increased because of the increased center-of- 
mass (CMS) energy y/s of the process. Due to the second point, the number of parameters 
on which the production rate depends is reduced since only the couplings of iV's to charged 
currents contribute (fiWN and eWN couplings). The second point does not apply in the 
case of fi~ fi + colliders, the case which we therefore do not consider in the present paper. 

As in the previous paper ||, our analysis is rather general, in the sense that we do not 
restrict ourselves to any specific (classes of) models. Further, we use basically the same set 
of programs developed by us there, this time adjusted to the case of efi colliders. We stress 
that, in contrast to the available literature, our numerical calculations include the effects 
of off-shell (nOS) intermediate iV's on the cross sections.[] This enables us to investigate 
deviations from the previously known cross sections (with two on-shell intermediate iV's), 
in the "20S" kinematic region (y/s > 1M > 2M\y) where both intermediate iV's can, but 
need not, be on shell (OS) - these deviations are often referred to as "finite width effects". 
Further, our general expressions allow us to calculate the cross sections also in the "10S" 



1 This was the case also in our previous paper ||. 



kinematic region (2M > ^/s> M + My/) where at most one intermediate N can (but need 
not) be on shell, and even in the "nOS" region (M+Mw > \/s) where both intermediate iV's 
always have to be off shell (nOS). For various choices of physical and kinematic parameters, 
we calculate the total cross section for the mentioned lepton-number violating reactions and 
the corresponding expected numbers of events for e/i colliders. 

II. REACTION AMPLITUDE AND ITS SQUARE 

Our starting point is the rather general Lagrangian density for the couplings of the heavy 
Majorana neutrinos N with W's and light leptons £j (j = 1, 2, 3; £~[ = e~, £^~ = /i~, £ 3 = r~): 

CnM*) =-AE BJplfrWil - jb)N{x)W^ + h.c. , (1) 

where B[ s can be regarded, at first, as basically free parameters (later we discuss some 
possible constraints on B[ s). In Eq. ([l]), g is the standard SU{2)l gauge coupling pa- 
rameter. We mention that the N£W couplings could include, in principle, the right-handed 
parts, but this possibility will not be considered here. Our choice would thus appear to 
suggest that the considered N (mass eigenstate) is primarily sequential, i.e., with the stan- 
dard SU{2) L x U(l)y assignments. However, our choice can approximately describe also 
many other scenarios, including the one in which N has an admixture of SU (2) i-sing\et 
component, provided the right-handed coupling analogous to (|T|) is reasonably suppressed. 
Further, it turns out that parameters B[ s affect the final results only via the combinations 

H'=\Bf*Bf\, H = J2\ B l ] \ 2 - ( 2 ) 

3=1 

We will restrict ourselves to considering the reactions e~fi + — > NN — > W ± £fW ± £j(—+ 
jets + £f£j), which represent an explicit violation of the total lepton number - characteristic 
for the intermediate heavy Majorana neutrinos. The reaction involves only the t+u (shortly: 
tu) channel - cf. Fig. [l|. The s channel (Z-mediated) is not allowed. 

For the calculation of the invariant amplitude M.f% (shortly: M) for the tu channel^ we 
used the 4-component spinors u^ a \q) = u(qa) and v^ a \q)=v(qa) as defined by Itzykson and 
Zuber || (I.Z.), but with the normalization convention as given in ||, i.e., u = y/2rrifUi.z. 
and v— y/2irtfVi,z.i where m/ is the mass of the fermion. Further, for the spin component 
quantum numbers of spinors, we use the following notational convention for tildes: a = 
1, 2 •£» a = 2, 1. The tu-channel amplitude M.^ is 

, n AMA (tu "> 

[{p-pe-Pwr-M^+iT w M w \ 

{[u(pf, a e )(^ w + 2p r e)Y(l-l5)u{p; a)] pfeab^C 1 +75^(1+75 Mp^ 
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u(p e ; a e )£(l-j 5 )u(p; a)] [v(p; a)(f w f+2p r e)(l+^ 5 )v(p e ; a e ) }+■■■ (3) 



M fi is defined by g: (fin\S\in) = \(2^)H^(P in - P fin )M fi , for \fin)^ 
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Here we use the notations of Fig. [I]. Further, </[ = g^7^; £;, = e^(p w ) and £„ = £\?'{p w ) are 
the real polarization vectors || of the final W's, with polarizations A, A = 1, 2, 3. M is the 
mass of iV's. Parameter A^ in Eq. ([3]) contains the strengths of the couplings of the two 
Majorana neutrinos in the tu channel 

A™ = 9 -B^B^BfB^Xu , (4) 

where Am is the phase factor in the Fourier decomposition of the Majorana field N{x) (cf. J7|; 
|Am| 2 = 1)- Pn hi Eq. (0) is the (scalar) denominator of the propagator of iV 

P N (PiP w ) = l/[( P i+p w ) 2 -M 2 + iMT N ] , (5) 

where T^ is the total decay width of N. The dots in Eq. fl3|) at the end stand for three 
analogous terms, obtained from the expression explicitly written above by replacements: (I) 
(p w ,£) <->• (jp w , e); (II) (pe, a) <r* (jj e ,a) and overall factor (—1); (III) combined replacements 
(I) and (II). In expression ([5[) we neglected the masses of all the charged leptons appearing 
(e,/i,r). 

We can reexpress any of the above terms in M.( tu > in alternative forms, by applying the 
following general identities: 

-rfu(qaY = (-l) a v(qa), -rfv(qa)* = (-lfu(qa), (6) 

and (4P) T = — 7°7 2 W)7 2 7°, where the Dirac basis and the conventions of || are used for 
7 M 's. We use (U), for example, if we want to employ, in scalar expressions in square brackets 
of Eq. (|3|), u(p; a) and u(p e ; cfy) instead of v(p; a) and v(p e ; a). Such transformations are very 
convenient when calculating (|.M| 2 ) = (|A^^ U ^| 2 ). Using the afore-mentioned tranformations 
judiciously, we can always end up with traces of 4 x 4 matrices involving u(q, /3)u(q, (3) = <fi 
and/or v(q, /3)v(q, (3) = $ (q = p,p, . . .; (3 = a, a, . . . — 1, 2). We performed these traces by 
writing programs in Mathematica M and employing there subroutine Tracer 0. Single 
traces involve up to fourteen 7 M matrices and (1 — 75) matrix. 

The expressions for (|.M| 2 ), obtained in the way described above, were then rewritten 
with short notations for various scalar products and for the contracted Levi-Civita tensors. 
They were then translated into an optimized FORTRAN form by employing Maple [JTOj , and 
fed into our main program for calculation of the total cross sections a. The integrand (|7Vf | 2 ) 
is a long expression, extending over tens of pages when printed out. The main program uses 



subroutines rambo |TT| and vegas [T2[ for numerical integration 



III. NUMERICAL CALCULATIONS AND RESULTS 

We stress that our general program can now deal with all kinematic situations: 
1.) those where both intermediate iV's can (but need not) be on shell [20S kinematic region: 
y/s>2M>2M w ]\ 

2.) those where only at most one N can (but need not) be on shell [10S kinematic region: 
2M>^s>M+M w }; 
3.) those where both iV's must be off shell [nOS kinematic region: M J cMw>\fs\- 



However, by the same methods as described above, we also calculated the amplitudes 
M. and their averaged squares (|.M| 2 ) when one N, or both TV's, are explicitly put on 
shell (10S, 20S expressions, respectively). These expressions turn out to be much shorter 
than the general (nOS) expressions with both iV's propagating. Analogously as for the 
nOS expressions (cf. previous Section), also the 10S and 20S expressions for (|.M| 2 ) were 
calculated with help of the subroutine Tracer and were then fed into separate FORTRAN 
programs where they were combined with the rambo and vegas subroutines. These (10S, 
20S) programs allow us to calculate the sum of cross sections for the considered reactions 
e~fi + — ► NN — > W + £~W + £~ (i,j = 1,2,3) where now one, or both, intermediate iV's 
are explicitly on shell. In fact, the 10S expression for (|A^| 2 ), for the sum of reactions 
e~/i + — > NN(OS) — > W + £~ N(OS) (2 = 1,2, 3), is multiplied simply by the branching ratio 
Br for the sum of the decay modes N(OS) — > W + £j (j = 1,2,3); the 20S expression for 
the reaction e"/i+ -> N(OS)N(OS) is multiplied by (Br) 2 . 

For numerical calculations, the input were values of y/s, M, and of the "mixing" parame- 
ters H' and H [cf. Eqs. (0)-(@)]- Parameter H' is a combined measure of the strengths of the 
eWN and fj,WN couplings and it affects the (tu) amplitude crucially (oc H'). Parameter H 
affects the total (|A^| 2 ) which is then formally proportional to H 2 (once H' is fixed). This 
is so because we sum over all the flavors of the two final light leptons £~ , £~ . All in all, 
in (|-M| 2 ) we thus average over the spin components (a, a) of the initial e~ and /i + (factor 
1/4), sum over the polarizations (A, A) of the two final PF + 's and over the spin components 
(ae,ai) and the flavors (i,j = 1,2,3) of the two final light leptons. In the general (nOS) 
expression, we have to include an additional factor 1/2 to avoid double-counting of the two 
W + 's, and factor 1/2 to avoid double-counting of the final leptons (i.e., of £^£^, or of twice 
£~£j when i ^ j). The (|A^| 2 ) thus defined can then be integrated over the entire phase 
space of the final particles (massive W's and massless £), and there is no double- count ingQ 
in the obtained total o. 

Further, the total decay width Y ^ of iV's, appearing in the denominators of their prop- 
agators, was determined by us at the tree level and under the assumption that the only 
(dominant) decay modes are N — > W £j (j = 1,2,3). Incidentally, this means that the 
branching ratio for the sum of the decay modes N — > W + £j (j = 1, 2, 3) is simply Br = 1/2, 
the value that is then also used in the 10S and 20S programs for cross sections (cf. discus- 
sion above). In this framework, we have Tn oc H, i.e., H is the decay width parameter of 
the iV's. 

In large classes of models, in which heavy neutrinos are sequential or have exotic SU(2)x 
U(l) assignments, the values of the "mixing" parameters H' and H [cf. Eq. (0)] are severely 
restricted by available experimental data (LEP and low-energy data) [TBl - fH] ]: 



3 In the simpler 10S and 20S expressions for (|.M| 2 ), these factors are different. For the 20S 
expression, the factor is (1/2) because the two intermediate on-shell iV's are indistinguishable, 
and thus their decay product pairs (W + £~ , W + £~) are also kinematically indistinguishable. For 
the 10S expression, in the 10S kinematic region, all four final particles are distinguishable via 
kinematics (i.e., invariant masses) of W + £~ pairs, and the factor is 1; in the 20S kinematic region, 
the two decay product pairs (W + £~ , W + £j) become kinematically indistinguishable, and the factor 
is (1/2). 



H' £ 0.0155 , H « 0.122 . (7) 

For orientation, we used for the values of H' and H these two upper bounds. The numerical 
results for the cross sections a of the reactions e~fi + — > NN — > W + £~W + £~ (summed over 
i, j = 1, 2, 3) are depicted in Figs. 0-|6|. 

In Figs. [||], the cross sections a are presented as functions of the Majorana neutrino mass 
M, at five different CMS reaction energies \fs = 0.5, 1, 2, 6 TeV, respectively. Calculations 
were performed in various kinematic regions (nOS, 10S, 20S regions) with the general 
(nOS) program which has the propagating intermediate iV's with finite width [cf. Eqs. ([3|), 
(I)]. Further, in the 10S and 20S kinematic regions, the cross sections were calculated also 
with the 10S program (in which one N is explicitly on shell). In the 20S kinematic region, 
the cross sections were performed in addition with the 20S program (in which both iV's are 
explicitly on shell). Comparing the results of the general nOS program with those of the 10S 
and 20S programs, we see from these Figures that the contributions of the configurations 
with off-shell Majorana neutrinos become important only at high ^Js > 2 TeV. For example, 
when ^/s — 2.0 TeV, the results of the nOS program in the 20S kinematic region (M < y/s/2) 
differ from those of the 20S (and 10S) programs already by up to 30%. One of the reasons 
for this trend lies in the fact that higher y/s (at a given M) mean that more off-shell 
configurations are available for the intermediate iV's. Further, in Figs. |2|4j] we see the slope 
increase associated with the change of kinematic regions: at M m y/s/2 (onset of the 20S 
region). In Fig. [5], only the 20S kinematic region is included. 

In connection with Figs. |2|-[5], one technical difficulty with the general (nOS) program 
should be mentioned. The results of the nOS program were not depicted in the lowest parts 
of the 20S kinematic regions (M much lower than yfs/2). In these regions of low M, the nOS 
program requires exceedingly high statistics to avoid large numerical uncertainties. These 
numerical uncertainties arise because, for low M's (~M^) only a very limited part of the 
rambo-generated phase space of final particles contributes to the cross section (that would 
be true even when T^ were, artificially by hand, prevented from approaching the zero value 
when M ^M\y). An additional important reason for these uncertainties lies in the fact that 
the integrand becomes very singularly peaked in small subregions of that mentioned part of 
the phase space. This is so because V ^ becomes small for small M's, and consequently the 
absolute square of expression @ is then almost equal to [k / (MT n)]5[{ph+Pw) 2 — M 2 ]. 

In Fig. ^| we plot the cross sections as function of the CSM energy y/s, at fixed chosen 
values of the Majorana mass M = 200, 400 and 600 GeV. Here we see that, at a given 
Majorana neutrino mass M, the full cross section usually (but not always) increases with 
increasing y/s. At high CSM energies y/s>2M (20S kinematic region), the cross section 
changes only slowly with increasing y/s. In Fig. |6|, for M = 200 GeV only results of the 20S 
program (in which both iV's are put on shell) are displayed, since the general (nOS) and 
the 10S programs show up the afore-mentioned numerical instabilities and large numerical 
uncertainties at so low values of M. Even for M = 400 GeV, the nOS program has significant 
numerical uncertainties (~10%) for a in the region y/s = 1-2 TeV, but these results were 
included nonetheless. 

One may ask how the depicted results change once we change the values of the "mixing" 
parameters H' and H [cf. Eqs. (|2|), (|7|)]. It turns out that, for the 20S kinematic region 
(y/s > 2M > 2M W ), the cross sections (and distributions) as given by the general (nOS) 
program depend on the parameter H only weakly. This can be understood first in the 

6 



formal way, because the absolute square of each of the two intermediate iV-propagators 
(when Tn^M) is approximately proportional to T^ociJ -1 (cf. previous discussion), while 
the sum of the squares of the amplitudes without iV-propagators is proportional to H 2 - 
the latter because of the two decay vertices N — > W + £~ and N — >■ W + £^ (cf. Fig. [I]) 
which are squared and summed (i,j = 1,2,3). This approximate if-independence can be 
expected on physical grounds: These cross sections are, for not too high -y/i, dominated 
by the configurations with two on-shell iV's; but the absolute squares of amplitudes for 
e~fi + —* N(OS)N(OS) [combined with the decay branching ratios (1/2) 2 for the N(OSys] 
are exactly if -independent. Parameter H, being responsible for the finite width (note: 
r^r oc H), is apparently one of the two principal parameters responsible in the 20S kinematic 
region for the deviation of the cross section from the pure 20S cross section (the other is 
a/s). On the other hand, in the 10S kinematic region (2M > ^/s > M + Mw), the in- 
dependence of the cross section becomes quite strong (approximately ocii), and in the nOS 
region (M+Mw>\fs) even more so (ocii 2 ). 

On the other hand, the dependence of the cross sections on the parameter H' [cf. Eqs. (|J), 
(||)] is always strong (ocii' 2 ). 

IV. SUMMARY AND CONCLUSIONS 

We investigated the reactions e~/i + — > NN — > W + £~W + £'~ where the intermediate iV's 
are heavy Majorana neutrinos (with mass M~10 2 GeV) and £,£' are light charged leptons 
(e,fj,,T). These processes violate the total lepton number and, if detected, would represent 
a clear signal of the Majorana character of (heavy) neutrinos. In contrast to the expressions 
available in the literature so far where the intermediate iV's are taken to be on shell, we 
calculated the amplitudes of these processes by employing finite width propagators of the 
iV's. This enabled us to calculate numerically the full cross section for the reactions, i.e., 
including the effects of the off-shell intermediate iV's. In contrast to the analogous reactions 
in e~e + (and \f~ /i + ) colliders, the considered reactions do not include a (Z-mediated) s 
channel, but only a (P^-mediated) tu channel, thus allowing us to reduce the number of 
relevant unknown coupling parameters. If the CMS process energies i/s are large (> 2 
TeV), the effects of the off-shell intermediate iV's can be significant even in the 20S kinematic 
region, i.e., in the region \/s> 1M (> 2Myy) where the configurations with both N's being 
on shell are allowed. Further, the number of such events at future e/x colliders would be 
in general high enough for detection. For example, if assuming the integrated luminosity 
10 4 pb" 1 , and Majorana mass M = 200 GeV (400 GeV), we would get about 10, 24, 29, 31 
(0, 27, 32, 188) events at ^=0.5, 1.0, 2.0 and 6.0 TeV, respectively. These numbers have to 
be multiplied by two if we include also the events with two positively charged leptons in the 
final state (the latter events also violate the total lepton number conservation, of course). 
We should stress that we used the highly suppressed values of the "mixing" parameters 
H' = 0.0155 and H = 0.122. These values are the maximal values allowed in theories where 
the heavy neutrinos are either sequential or have exotic SU(2) x U(l) assignments. If these 
bounds do not apply, the number of events can increase dramatically since it is proportional 
to H' 2 . 

In the paper, we ignored the questions connected with the experimental difficulties of 
detecting the discussed process unambiguously. In particular, there are problems connected 



with the identification of the (on-shell) W's and some of the light charged leptons (t's). 
Further, we ignored the possibility that the Majorana neutrino mass M is very low - below 
M\y (the jy + 's are then off shell) |TI|. In such interesting cases of low M, however, additional 



experimental problems would arise in the identification of the process since the two W /+ 's 
are now intermediate off-shell particles. 
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FIG. 1. Reaction e~(p, a)/i + (p, a) -> N*N* -> W+ (p w , X)W+ {p w , X)£^ {p £l at)£~ (p e ,a e ). Attaching 
the legs of the four final particles to the two -/V-propagators differently leads to the other contributing 
(tw-channcl) diagrams. The CMS squared energy of the process is s = (p+J5) 2 . 
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FIG. 2. Sum of the total cross sections for the reactions e~^+ — > NN — > W + W + £^iJ (i,j — 1,2,3; 
£i=e, £2 = ^, ^3 = T ), as function of the Majorana neutrino mass M, for the CMS process energy ^=0.5 
TeV. Summation over the final, and average over the initial polarizations were carried out. The values of 
the "mixing" parameters H' and H of Eq. (0) were set equal to the upper bounds (R). The results of the 
general (nOS) program are displayed as full lines, in various kinematic regions. The results of the simpler 
10S program are displayed as triangles, in the 10S and partly in the 20S kinematic regions. The results of 
the short 20S program are displayed as crosses, in the 20S kinematic region. 
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FIG. 3. Same as Fig. |, but for ^=1.0 TcV. 
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FIG. 5. Same as Fig. |, but for ^=6.0 TeV. 
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FIG. 6. Sum of the mentioned total cross sections for the reactions 



i+ -► NN -► W+W + £-£: 



(i,j = 1,2,3; £i = e, £2 = ^, ^3 = T ), as function of the CMS process energy y/s, for Majorana neutrino mass 
M = 600, 400 and 200 GeV. 
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